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EXPERIMENTAL AND CALCULATED SUPERSONIC 
FLUTTER CHARACTERISTICS OF MODELS 
OF THE X-15 HORIZONTAL AND 
VERTICAL TAIIS* 

By William T. Lauten, Jr., and Robert W. Hess 


SUMMARY I 


A series of tests have been conducted at Mach numbers of 1.64, 2.0, 
3 . 0 , and 4.0 to investigate the possibility of flutter on dynamically 
and elastically scaled models of the all-movable horizontal- and vertical- 
tail surfaces of the X-15 airplane. No flutter was obtained on the scaled 
models . 


Flutter was obtained on four weakened models of the horizontal tail 
at Mach numbers of 1.64, 2,0, and 3*0 and one weakened vertical -tail model 
at a Mach number of 3-0* The test results and structural information are 
presented. Also presented are the results of flutter calculations using 
aerodynamic forces derived from piston theory. The agreement between cal- 
culated and experimental flutter speeds is very good at a Mach number 
of 3.0 for the horizontal tail, but at the lower Mach numbers the 
ment is less satisfactory. 


INTRODUCTION 


agree-. . 


Flutter characteristics of various components of the X-15 airplane 
have been determined at the Langley Research Center. As part of this 
program, tests have been made on dynamically and elastically scaled 
models of the X-15 components in the various speed ranges of the airplane. 
Previous tests on the X-15 components for hypersonic speeds are reported 
in references 1 and 2 . 




This paper presents the results of tests in the Langley 9 - by l 8 -inch 
supersonic aeroelasticity tunnel of models of the horizontal- and vertical- 
tail surfaces at Mach numbers of 1.3, 2.0, 3.0, and 3 . 98 . Since none of 
the scale models fluttered, the results that are presented are for model 
surfaces whose stiffness has been reduced until flutter occurred. A com- 
parison is made with calculations in which the aerodynamic forces are 
derived from piston theory. (See ref. 3*) The first three natural modes 
of the horizontal -tail models and the first two natural modes of the 
vertical-tail model were used in the calculation and were experimentally 
determined by the method outlined in reference 4. 
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SYMBOLS 


speed of sound, ft /sec 
root semichord of model, ft 
local semichord, ft 
frequency, cps 

mode shape, Displacement of any point on surface 

’ Displacement of point of maximum displacement 

span of model, ft 

Mach number 

mass of model segments, slugs 
mass of increments of model, slugs 


model to air mass density ratio, 


frequency, radians / sec 


Mass of wing 


r 

Jo 


jtp / c^dy 


dynamic pressure, lb/sq ft 




P 


air density, slugs/cu ft 
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V velocity, fps 

y spanwise coordinate 


Subscripts: 


1 , 2 , 3, 4 

e 

c 


indicate natural modes in ascending order of frequency 
experimental flutter data 
calculated flutter data 


MODEL DESCRIPTION 


All models were supplied by the manufacturer. Although the models 
were of similar construction, they were scaled in stiffness and mass on 
the basis of the dynamic pressure available in the test facility. Since 
the available dynamic pressure decreases with an increase in Mach number 
above M = 2, the models designed for the higher Mach numbers were lighter 
and more flexible as the Mach number increased above M = 2. 

The models are identified in the text by the following scheme: HT 

and VT designate horizontal- and vertical-tail, respectively. The 
suffix contains the model number and the letter H or C. The stiff ix C 
designates those models which represent the airplane in the cold condi- 
tion, whereas the suffix H designates the models which have been weakened 
to represent the airplane in the hot condition. For example, a model 
designated VT-2 is vertical-tail model 2 whereas model HT-4H is horizontal- 
tail model 4 which has been weakened to represent the airplane in the hot 
condition. 


Horizontal Tail 

The models of the horizontal stabilizer were l/l2 scale with an 
exposed surface aspect ratio of approximately 2.5, a taper ratio of 0.305, 
and a sweep angle of 45° at the quarter-chord line. The airfoil section 
was an NACA 66AOO5 modified so that it was 1 percent thick at the trailing 
edge with a straight-line fairing to the point of tangency. The airfoil 
ordinates are listed in table I. 

Of the five horizontal-tail models of this investigation, four were 
designed to simulate surfaces reduced in stiffness by aerodynamic heating. 
The stiffness of the fifth model was not reduced and represented the case 
of the airplane not weakened by aerodynamic heating. Since no flutter 
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was obtained on the dynamically |cal^J. qip^els^ the spindle spring 
restraint was reduced below the scaled values until flutter occurred. 

All the horizontal-tail models were built with an aluminum box spar 
and rib construction. Small lead weights were distributed over the plan 
form in order to obtain the proper dynamic characteristics. To obtain 
the airfoil shape the section was filled with balsa and covered with 
lightweight silk. Figure 1 is an X-ray picture which shows typical 
construction. 

A top-view drawing of the model mounted in the test section is 
shown in figure 2( a) . The model base block served as a model mount and 
also as a spacer to support the model in the airstream beyond the tunnel 
boundary layer. A reflection plane, also shown in figure 2(a), was 
attached to the model base block just inboard of the wing root. The 
model was supported in its base block by means of two flexure pivot 
springs which simulated the spindle bearings of the prototype. The pivot 
was at 25 percent of the mean aerodynamic chord. A spring was mounted 
ahead of the model spindle and attached to the spindle with a screw as 
is shown in figure 2(b). The stiffness of the pitch spring and the pivot 
support combined to simulate the actuator stiffness of the prototype. 


Vertical Tail 

0 

The vertical- tail models were l/l2 scale with an aspect ratio 
of 0.359 and- a leading-edge sweep angle of 29°. The airfoil section was 
a wedge with an enclosed angle of 10°. Figure 3 shows top- and end-view 
drawings of the model as mounted in the test section. The vertical tail 
was built in two parts . Approximately one-fourth of the span was rigidly 
attached to the base block. On one model speed brakes comprised some 
34 percent of the rear portion of the inboard section. Figure 4 is a 
photograph of the trailing edge of this model showing the speed brakes 
extended. The movable portion of the tail was attached to a spar which 
was attached to both the fixed portion and the model base block with 
springs . 

Two models of the vertical tail were tested; one constructed with 
a magnesium frame covered with sheet balsa and silk did not flutter. 

This model had the speed brakes and is shown in figure 4. The second 
model which did flutter was much lighter and weaker. An X-ray photograph 
of this model is shown in figure 5- The spax was formed by wrapping a 
balsa box with aluminum-alloy sheet. The tail was then attached to the 
supporting spindle at the root-spar cap which was integral with the model 
and may be seen as the daxk area at the root of the movable portion of # 

the vertical tail. The fixed portion of the second model was not equipped 
with speed brakes. 
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INSTRUMENTATION 


y 


A recording oscillograph was used to obtain a continuous record 
during each test of tunnel stagnation, temperature, and pressure. Signals 
from strain gages mounted on the model spar were recorded on the same 
record and were used to determine flutter onset and frequency. Motion 
pictures, at 1,000 frames per second, were taken of the model during each 
tunnel run. 


LABORATORY MEASUREMENTS 


Prior to each run the frequencies and node lines for the first four 
natural modes were determined for the horizontal tail and the first two 
modes for the vertical tail. Figures 6(a) and 6(b) show typical sets of 
node lines. 

The mode shapes were determined by the method of reference k on at 
least one, usually the stiffest wing in each series of tests. These 
mode shapes are presented in tables II to VII. 


Horizontal Tail 

After flutter testing was completed, the horizontal-tail models which 
were damaged but not completely destroyed were cut up in order to deter- 
mine the section properties. The pitch inertia, roll inertia (about the 
strip center of gravity), and mass of the strips shown in figure 7 are 
given in table VIII. The section mass moment of inertia was determined 
by use of a bifilar suspension. The mass distribution in each section 
was also determined. Since the small lead weights cause discontinuities 
in the chordwise mass distribution, the streamwise sections were cut into 
small blocks. The apparently random distribution of cuts shown in fig- 
ure 7 was decided upon to minimize the amount of mass lost in cutting. 
Table IX lists the mass of the blocks shown in figure 7 but does not 
include the mass of material lost in cutting. 

Although the model designated as HT-7C was destroyed during the 
testing, its construction was so nearly identical to that of model HT-1H 
that a mass distribution for HT-7C was generated by consideration of the 
differences in the location and mass of the lead weights in the two models 
and the total mass of the model determined before testing. This mass 
distribution is also given in table IX. 
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Vertical Tail 

It was necessary to calculate the mass distribution of the vertical- 
tail model VT-2 since it was destroyed by flutter. The thickness of the 
aluminum-alloy- spar skin, 0 .005 inch, and the weight of the lead ballast 
and spar attachment cap were known. (See fig. 5*) The balsa sheet was 
assumed to be 0.0625 inch thick with a density of 6 pounds per cubic foot. 
The total weight of the vertical-tail model as determined from the summa- 
tion of the measured and calculated weights of the various parts was within 
2 percent of the total measured weight of the model. 

Chordwise mass distribution curves, based on the previously mentioned 
assumptions, are given in figure 8; The ordinates of these curves are in 
slugs per square foot of control surface erea and axe based on the thick- 
ness and density of material at a particular point. For 10-percent span 
to 90-P e ncent span the values for that portion of the chord which includes 
the spar have been adjusted to include the mass of the spar webs. The 
mass of the trailing edge, the mass of the lead balance weights, and the 
mass of the spar cap are not included. Neither is the mass of the two 
intermediate ribs included. The masses of the balance weights and the 
spar caps are shown in figure 5- The mass of increments of the two 
intermediate ribs and the trailing edge is given in table X. 

Since these discrete masses are difficult to represent on a mass 
distribution curve, they are added into the flutter calculations as an 
incremental correction after the generalized mass of the model is cal- 
culated by utilizing the mass distribution curves. 


COMPARISON OF FREQUENCY SPECTRUM, MODE SHAPES, AND 
NODAL PATTERNS FOR MODELS AND PROTOTYPE 


Information on the tail surfaces of the prototype has become avail- 
able and is compared with the vibration characteristics of the models. 
The full-scale information was obtained in vibration surveys made pre- 
liminary to some tests of full-scale tail surfaces in the Langley by 
6-foot thermal structures tunnel. 

A comparison of the nodal patterns and frequency spectrums of the 
models and prototype are shown in figure 6(a) for the horizontal tail 
and in figure 6(b). for the vertical tail. A graphical comparison of the 
second mode of models HT-7C and HT-8H with the second mode of the proto- 
type is shown in figure 9. The prototype modes which were measured are 
presented in tabular form in table XI. 
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TEST RESULTS 


A summary of the test results is listed in table XII. The table 
lists the tunnel conditions at flutter. For model HT-8H, where no flutter 
was obtained on the weakest feasible configuration, maximum q is listed. 
Also listed are the density ratio of model mass to air mass |i, the model 
natural frequencies, the flutter frequency, and the stiffness-altitude 
bav> 

parameter \[jl, where a >2 is the frequency of the second natural mode 

£L 

in radians per second. The results of the flutter tests are also pre- 
sented graphically in figures 10(a) and 10(b) which are plots of the 
stiffness-altitude parameter against Mach number for the horizontal and 
vertical tails, respectively. In these figures a constant dynamic pres- 
sure is a straight line through the origin. The flutter region is below 
and to the right of the flutter boundary, larger values of the stiffness- 
altitude parameter indicating less likelihood of flutter. Also shown is 
an airplane flight line based on a constant dynamic pressure of 2,500 pounds 
per square foot. 


Horizontal Tail 

No flutter was obtained on any of the dynamically scaled models of 
the horizontal tail at the Mach number for which they are designed. In 
order to obtain flutter information, the spring restraint was reduced 
below scaled values until flutter occurred. A comparison of the airplane 
flight path and the experimental results, plotted in figure 10(a), indi- 
cates that the flutter margin established by the models is smallest at 
M = 1.64. Since the stiffness-altitude parameter is inversely propor- 
tional to the square root of the dynamic pressure, it can be noted that 
at this Mach number the dynamic pressure corresponding to the experimental 
flutter point is 60 percent greater than the dynamic pressure corresponding 
to the flight path. 


Vertical Tail 

No flutter was obtained on the first vertical tail tested, model VT-1. 
The tests were partially inconclusive in that the speed brakes of the model 
were bottomed out on the stops by the aerodynamic loading before the 
desired dynamic pressure was obtained. This result was due to the fact 
that the spring stiffness of the speed brake was determined from frequency 
considerations and that the spring was not large enough to be preloaded 
to a point where it could withstand the aerodynamic loads. The tests 
were extended to the maximum available dynamic pressure with the weakest 
feasible root restraint in an effort to obtain flutter on the all-movable 
portion of the surface. 
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After no f^utt^r was ol|tjiijied on the first model, a second model VT-2 
was constructed, as described in the section "Model Description," and 
tested at M = 3. This model fluttered to destruction; the results are 
presented graphically in figure 10(b). At a Mach number of 3 the experi- 

bCJDp r— 

mental flutter point is at a — - Jn of 3*5 whereas the flight boundary 

EL 1 

b(no i— 

is at a \m of 9; thus, a very large margin of safety is indicated. 

EL * 


CALCULATED RESULTS 


Modal-type flutter analyses were made by using the method of cal- 
culation set forth in reference 5* Aerodynamic forces were based on 
piston theory (ref. 3) which assumes that the local pressure on one point 
of the wing is a function only of the downwash velocity at that point and 
is independent of disturbances at other points of the wing. The first 
three natural modes of the horizontal tail and the first two natural modes 
of the vertical tail were used in the calculations, the mode shapes having 
been determined by the method of reference 4. The results of these cal- 
culations are summarized in table XI and are compared graphically with 
the experimental results in figures 10(a) and 10(b). 


Horizontal Tail 

The agreement between experimental and calculated flutter speeds 
for the horizontal tail at M = 3*0 was very good, with a difference 
of only 1 percent in the conservative direction. As might be expected 
when using piston theory, the agreement at the lower Mach numbers was 
not as good. The difference between the two values was 11 percent at 
M = 2.0 and l6 percent at M = 1.6^, both unconservative. Although 
no flutter was obtained at M = 3*98> in view of the agreement between 
experimental a n d calculated values at M = 3*0, it was felt to be worth- 
while to make similar calculations for the model tested at M = 3*98* 

These calculations showed that flutter should occur at a higher dynamic 
pressure than was available at this Mach number. 

As a matter of interest, and as a check on the effects of possible 
inaccuracies in the determination of mode shape, the mode shape of the 
second mode of the horizontal-tail model, the mode to which these partic- 
ular calculations are most sensitive, was arbitrarily changed by decreasing 
the deflection of the trailing edge by 20 percent and the deflection of 
the 75-P e rcent chord line 10 percent. This change in mode shape is felt 
to be greater than any expected experimental inaccuracies. The resultant 
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calculated flutter speed at a Mach number of 3-0 was 3.4 percent lower 
than the original calculated value. 


A second alternate calculation was made by using a less exact mass 
distribution. Instead of several mass increments along each streamwise 
strip, only three were used. These were the portion of the strip ahead 
of the spar, the spar section, and the portion of the strip behind the 
spar. The mode shapes used were the same as those for the original cal- 
culation. This modification resulted in a calculated speed which was 
6.3 percent lower than the original value. 

A third alternate calculation was made by omitting the third mode 
from the calculations. This variation in the input data resulted in a 
calculated speed which was 7*5 percent higher than the original value. 

It is interesting to note that for this particular model the intro- 
duction of the first two variations (change in mode shape and mass dis- 
tribution), which might occur as a result of experimental errors or 
approximation, result in a more conservative flutter speed. The omission 
of the higher frequency mode raises the calculated flutter speed and thus 
tends to be unconservative. 


Vertical Tail 

The results of the flutter calculations for model VT-2 are presented 
in table XII and figure 10(b). It may be seen that the calculated values 
are 28 percent conservative. There is no obvious reason why the differ- 
ence between experimental and calculated results for the vertical tail is 
so much larger than the difference between the same two quantities for 
the horizontal tail. The lower aspect ratio and the unconventional air- 
foil shape of this configuration may be contributing factors. 


CONCLUSIONS 


Tests have been made on dynamically and elastically scaled models 
of the all-movable horizontal and vertical surfaces of the X-15 airplane 
to investigate the possibility of flutter. The models with design stiff- 
ness did not flutter. Flutter points were obtained at Mach numbers 
of 1.64, 2.0, and 3-0 on weakened models of the horizontal tail and at 
a Mach number of 3-0 on a weakened model of the vertical tail. 

1. The minimum difference in dynamic pressure between the tentative 
flight path and the experimental model flutter boundary of the two sur- 
faces occurs at a Mach number of 1.64. Based on the stiffness-altitude 
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parameter the flutter dynamic pressure is 60 percent greater than the 
expected dynamic pressure at this Mach number. 

2. The agreement between calculated and experimental flutter speeds 
was excellent for the horizontal tail at a Mach number of 3> but the cal- 
culated values tended to deviate, and become unconservative, from the 
experimental values at Mach numbers of 2 and 1.64. 

3. Flutter was only obtained on one model of the vertical tail at 
a Mach number of 3* This result indicated that the airplane vertical 
tail has a very large margin of safety from flutter at this Mach number. 

A comparison of the experimental with the calculated flutter speed showed 
that the calculated speed was 28 percent lower than the measured speed. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., August 19, 1959. 
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TABLE I.- ORDINATES FOR NACA 66A005 (MODIFIED) AIRFOIL 


[0 


Ordinates in percent of chord 


0 


y u 



X 

II 

0 

0 

.10 

.269 

.25 

.408 

.50 

• 531 

• 75 

.590 

1.25 

.650 

2.50 

• 791 

5.00 

1.048 

7.50 

1.270 

10.00 

1.460 

15.00 

1.766 

20.00 

2.001 

25.00 

2.182 

50.00 

2.518 

35.oo 

2.4 i 6 

40.00 

2.476 

45.00 

2.500 

50.00 

2.488 

55.00 

2.458 

60.00 

2.546 

65.OO 

2.176 

*67.00 

2.085 

100.00 

.500 


Straight-line fairing from 67- to 100-percent chord. 



TABLE II.- MODE SHAPE FOR MODEL HT-1H 
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TABLE III . - MODE SHAPE FOR MODEL HT-2H 
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TABLE V.- MODE SHAPE FOR MODEL HT-7C 
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TABLE VII.- MODE SHAPES FOR MODEL VT-2 
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TABUS IX.- MASS OF MODEL SBOMESTS SHOWN EN FIGURE 7 




AT TRAILING EDGE AND RIBS FOR VT-2 



l Access hole. 




































TABLE XI.- MODE SHAPES FOR X-15 PROTOTYPE HORIZONTAL TAIL 
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TABLE XII.- SUMMARY OF TEST RESULTS 



No flutter 










Figure 1.- X-ray photograph of model of horizontal tail showing typical construction 



Reflection plane 



(b) Photograph showing detail of model restraint. L-57-4258.1 








L-280 




photograph of model VT-2. L-59-602^ 





(a) Horizontal tail. 

Figure 6.- Comparison of nodal patterns and frequency spectrums of 
model and prototype tail surfaces. 
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Figure 8.- Cho: 
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Figure 10.- Plot of stiffness-altitude parameter against Mach number showing experimental 
flutter points, calculated results, and 2 , 500 q flight boundary. 








